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This  report  covers  the  period,  October  1,  1989  to  September  30,  1992.  It 
was  preceded  by  contract  number  N00014-84-K-0640  for  which  we  previously 
reported  completion  of  eight  journal  articles.  These  papers  were  mainly 
concerned  with  the  development  of  an  ice  model  and  the  coupling  of  that  model 
to  the  (now  standard)  Princeton  ocean  model.  The  studies  were  largely  based  on 
either  one-dimensional,  time  dependent  models  or  two-dimensional  (x,z) 
models.  The  latter  included  a  model  study  of  atmospheric-ice  interactions. 

The  more  current  grant  produced  four  papers.  The  first,  “One  hundred 
years  of  Arctic  ice  cover  variations  as  similated  by  a  one-dimensional  ice-ocean 
model"  was  an  application  of  our  time  dependent,  one-dimensional  model  to 
climate  time  scale  variability.  The  second  paper  "An  equation  of  state  for 
numericla  models  of  oceans  and  estuaries"  was  a  modification  of  the  UNESCO 
equation  of  state  which  rendered  it  useful  for  numerical  ocean  models. 

The  next  paper,  "Modeling  deep  convection  in  the  Greenland  Sea" 
represented  the  emergence  of  our  fully  three-dimensional,  ice-ocean  coupled 
model.  Here  it  was  applied  to  the  process  study  of  deep  convection  and  bottom 
water  formation.  Finally,  the  paper  "Modeling  the  seasonal  variability  of  a 
coupled  Arctic  ice-ocean  system"  was  a  modeling  study  of  the  Arctic  Ocean  and 
the  Norwegian-Greenland-Barents  Seas  and  included  multi-year  calculations  and 
sensitivity  studies.  The  first  page  from  these  papers  are  attached  to  this  report. 


□  □ 


FINAL  TECHNICAL  REPORT 
Grant  No.:  N00014-90-J-1268 


Page  2 


It  should  be  noted  that  this  grant  and  its  predecessors  contributed  to 
general  improvement  of  the  Princeton  Ocean  Model  which,  through  an  electronic 
users  group,  is  now  in  use  throughout  the  world;  this  includes  Naveocear.o  and 
other  Navy  sponsored  research  projects. 
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One  Hundred  Years  of  Arctic  Ice  Cover  Variations  as  Simulated 
by  a  One-Dimensionai,  Ice-Ocean  Model 

S.  Hakkinen  and  G.  L.  Mellor 


Princeton  University,  Princeton,  Sew  Jersey 


A  one-dimensional  ice  oeean  model  consisting  of  a  second  moment,  turbulent  closure,  mixed  layer 
model  and  a  three-layer  snow-ice  model  has  been  applied  to  the  simulation  of  Arctic  ice  mass  and 
mixed  layer  properties.  The  results  for  the  climatological  seasonal  cycle  are  discussed  first  and  include 
the  salt  and  heat  balance  in  the  upper  ocean.  The  coupled  model  is  then  applied  to  Ihe  period 
1880-1985,  using  the  surface  air  temperature  fluctuations  from  Hansen  el  ai.  (1983)  and  from  Wig'ev 
et  al.  (1981).  The  analysis  of  the  simulated  large  variations  of  the  Arctic  ice  mass  dunng  this  period 
(with  similar  changes  in  the  mixed  (aver  salinity)  shows  that  the  variability  in  the  summer  melt 
deternunes  to  a  high  degree  the  variability  in  the  average  ice  thickness.  The  annual  oceanic  heat  flux 
from  the  deep  ocean  and  the  maximum  freezing  rate  and  associated  nearly  constant  minimum  surface 
salinity  flux  did  not  vary  significantly  interannually.  This  also  implies  that  the  oceanic  influence  on  the 
Arctic  ice  mass  is  minimal  for  the  range  of  atmospheric  variability  tested. 


t.  Introduction 

The  main  goal  of  this  paper  is  to  explore  the  behavior  of 
the  ice-upper  ocean  system  on  climatic  time  scales  and 
especially  the  asymmetry  of  the  freezing  and  melting  sea¬ 
sons.  We  study  here  only  the  Arctic  ice  cover,  which  is  a 
prominent  climate  indicator  since  it  contains  a  large  multi¬ 
year  ice  component.  Also,  a  goal  is  to  further  test  the 
one-dimensional  ice-ocean  model  of  Mellor  and  Kantho 
[1939]  (hereinafter  referred  to  as  MK),  in  which  a  change  in 
a  boundary  condition  at  the  ice-ocean  interface  has  been 
implemented. 

The  influence  of  the  polar  ice  cover  on  climate  derives 
mainly  from  the  fact  that  the  ice  cover  can  insulate  the  ocean 
and  thereby  modify  the  surface  heat  flux  by  an  order  of 
magnitude.  The  insulating  effect  depends  on  the  reflective 
properties  of  the  ice  surface,  the  ice  concentration,  and  the 
ice  thickness  [Maykur,  1978,  1982;  Maykut  and  Unter- 
Steiner,  1971).  Climate  models  are  shown  to  have  a  strong 
sensitivity  to  snow  and  ice  cover  albedo  in  Polar  regions 
[Manabe  and  Stouffer,  1980). 

Leads  are  very  important  to  climate  response,  since  they 
release  heat  to  the  atmosphere  in  the  wintertime  and,  as  a 
secondary  effect,  lower  the  surface  albedo.  Heat  loss 
through  a  1%  lead  area  can  be  as  large  as  the  conductive  loss 
through  the  remaining  99%  of  the  area  covered  by  thick  ice. 
The  estimates  for  lead  area  range  from  1  to  5%  for  winter 
months.  Variability  in  the  lead  area  has  been  shown  to  affect 
considerably  the  surface.atr  temperatures  in  climate  models 
[Ledley,  1988 a,  b\. 

Monitoring  the  changes  in  ice  cover  have  been  limited  lo 
records  of  ice  extent  and  concentration  in  the  marginal  ice 
zones  before  the  advent  of  satellites  in  the  1970s  [Walsh  and 
Johnson,  1979)  or  to  a  local  ice  index  such  as  compiled  for 
the  Icelandic  coastal  waters  [ Sigtryggson ,  1972).  The  latter 
record  beginning  in  1900  shows  that  there  are  periods  with 
extremely  heavy  ice  around  1910  and  in  (he  end  of  the  1960s. 
Remote  sensing  techniques  have  now  made  it  possible  to 
monitor  the  ice  concentration  changes  from  space  (e.g.,  the 
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passive  microwave  sensors  such  as  electrically  scanning 
microwave  radiometer  (ESMR).  ( Parkinson  et  al.,  1987), 
scanning  multichannel  microwave  radiometer  (SMMR).  or 
special  sensor  microwave/imager  (SSM/I1).  The  ice  thick¬ 
ness  variations  pose  a  much  more  difficult  monitoring  task. 
Even  though  SMMR  and  SSM/I  can  resolve  thin  ice  and 
multiyear  ice  components  of  the  ice  concentration,  they 
C3nnot  measure  actual  thickness  of  the  sea  ice.  Upward 
looking  sonar  along  submarine  tracks  has  been  the  only 
measurement  able  to  give  a  large-scale  survey  of  the  ice 
thickness  variability  in  the  Arctic  [Bourke  and  Garrett,  1987; 
Wadhams,  1988). 

On  climatic  time  scales  the  ice  thickness  is  an  important 
factor  in  controlling  the  conductive  heat  flux  through  ice  and 
thus  the  salinity  flux  to  the  ocean  arising  from  bottom  ice 
accretion.  Because  ice  stores  and  releases  fresh  water,  the 
ice  thickness  and  concentration  variations,  which  determine 
the  ice  mass,  are  important  quantities  in  the  hydrological 
cycle.  (Note,  however  that  the  freshwater  reservoir  of  Arctic 
ice  is  a  thousand  times  smaller  than  that  of  the  Greenland  ice 
sheet  [Goody,  1980).  This  sea  ice-fresh  water  source  can 
affect  the  subpolar  gyre,  mainly  through  the  ice  exported 
from  the  Arctic  to  the  Greenland  Sea.  and  can  significantly 
modify  the  oceanic  surface  layers  in  the  areas  conducive  to 
water  mass  modification  sur  t  as  in  the  Greenland,  Iceland, 
and  Labrador  seas  [Swift  o>  d  Aagaard,  1981). 

The  ice  export  from  the  Arctic  through  the  Fram  Strait 
depends  also  on  the  winds  and  ocean  currents  in  the  outflow 
region.  Estimates  of  yearly  averaged  areal  outflows  are  0.9  x 
|04  km:  [Vowinckel,  1964],  0.84  x  106  kmJ  [Moritz.  1988), 
and  1.35  x  104  kmJ  ( Vinje  and  FinnekAsa,  1986J.  Large- 
scale  ice-ocean  models  have  predicted  the  ice  export  to  vary 
up  to  a  factor  of  3  interannually,  with  a  mean  of  1400  km J  for 
volumetric  and  a  mean  of  0.63  x  to4  km J  for  the  areal  export 
[ Walsh  et  al.,  1985).  The  ice  export  variability  in  the  model 
of  Walsh  et  at.  is  very  sensitive  to  the  wind  forcing,  which 
contributes  75%  to  the  ice  drift.  This  is  in  contrast  to 
estimates  where  wind  and  current  drift  portions  are  30%  and 
70%.  respectively,  as  derived  from  observed  buoy  drifts  and 
geostrophic  winds  [Moritz.  1988).  Empirical  models  for  the 
icc  outflow  based  on  the  air  pressure  difference  across  the 
Fram  Strait  suggest  that  the  wind-driven  flow  may  have 
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An  Equation  of  State  for  Numerical  Models  of  Oceans  and  Estuaries 

George  L.  Mellor 

Atmospheric  and  Oceanics  Sciences  Program.  Princeton  University.  Princeton.  Seer  Jersey 
26  October  1990  and  3  January  1991 


1.  Introduction 


Numerical  models  of  oceans  or  estuaries  require  an 
equation  of  state  in  order  to  relate  density,  salinity, 
temperature,  and  pressure.  The  needs  are  for  an  equa¬ 
tion  that  is  reasonably  efficient  numerically  and  has  a 
wide  range  of  application.  In  addition,  since  numerical 
models  solve  for  potential  temperature,  one  would  like 
an  equation  of  state  in  the  form 


p  =  p(5,  0,  p)  (1) 


where  p  is  density,  5  is  salinity,  0  is  potential  temper¬ 
ature,  and  p  is  pressure. 

In  the  following  we  will  denote  property  functions 
as  in  Eq.  ( 1 )  where,  for  example,  dp/dS  denotes  the 
partial  derivative  of  density  with  respect  to  5  holding 
S  and  p  constant.  We  denote  spatial  functions  with  a 
superposed  tilde  so  that,  for  example,  dp/dx  denotes 
the  partial  derivative  of  density  with  respect  to  x  hold¬ 
ing  y  and  z  constant. 

Models  require  horizontal  density  gradients,  (dp/ 
dx,  dp/dy),  which  are  readily^  obtained  from  Eq. 
(  1 )  and  known  spatial  fields  of  T,  S,  and  p.  They  also 
requirt  the  vertical  static  stability  obtained  either  from 


or  from 


g  p  dS  dz  dB  dz  _ 


g  p[dz  c,2  dz 


(2a) 


(2b) 


where  c,2  ■■  [(dp/dp),.*]-'  is  the  speed  of  sound 
squared.  The  term  jV2  governs  vertical  mixing  either 
through  a  Richardson  number  formulation  (e  g.,  Munk 
and  Anderson  1 948 )  or  by  a  turbulence  closure  scheme 
(e.g.,  Mellor  and  Yamada  1982).  Equations  (2a)  and 
(2b)  are  equivalent  since  dp/dz  -  (dp/dS)(dS/dz) 
+  (dp/d8)(d0/dz)  +  (dp/dp)(dp/dz).  However,  nu¬ 
merical  calculations  using  (2a)  or  (2b)  may  differ  as 
discussed  below. 


Corresponding  author  address:  Dr.  George  L.  Mellor,  Program  in 
Atmospheric  and  Oceanic  Sciences,  P.O,  Box  CN?  ■  J,  Sayre  Hall, 
Princeton  University,  Princeton,  NJ  08344 — 6571 . 


Bryan  and  Cox  (1972)  use  a  three-dimensional 
model  whose  vertical  coordinate  corresponds  to  con¬ 
stant  depth  levels.  To  obtain  density  they  preprocess  a 
table  of  coefficients  for  each  depth.  On  the  other  hand, 
we  have  a  three-dimensional,  primitive-equation  ocean 
model  that  uses  a  sigma  coordinate  system  in  the  ver¬ 
tical  (Oeyetal.  1986a.b;  Blumberg  and  Mellor  1987). 
It  has  been  applied  to  coastal  oceans  and  estuaries  but, 
increasingly,  is  also  applied  to  deep  water  ocean  do¬ 
mains.  Other  sigma  coordinate  models  have  appeared 
and  are  still  appearing  in  the  literature.  For  sigma  co¬ 
ordinate  models,  the  Bryan-Cox  scheme  is  inappro¬ 
priate. 


2.  Analysis 

The  now  standard,  UNESCO  (United  Nations  Ed¬ 
ucational,  Scientific,  and  Cultural  Organization) 
equation  of  state  is  fairly  expensive  computationally 
and  is  overly  precise  relative  to  the  capability  of  nu¬ 
merical  models  to  produce  comparable  precision.  Fur¬ 
thermore,  it  is  an  equation  of  the  form,  p  =  p(S,  T, 
p)\  i.e..  in  situ  temperature  is  an  independent  variable 
rather  than  potential  temperature. 

A  less  precise  but  substantial’.)'  simpler  formula  of 
the  desired  form  is 

p(S,e,p)  =  p<S.P,0)  +  p(l  -CpJ  (3) 

where  p( S,  0,  0)  can  be  taken  from  the  UNESCO  for¬ 
mula  ( the  complete  formula,  p(S,  T ,  p ),  increases  the 
computational  time  by  about  a  factor  of  3 )  since  ©  **  T 
for  p  =  0.  The  second  term  in  Eq.  (3),  the  pressure 
dependent  part,  contains  c  =  c(S,Q,p)  and  a  constant 
C.  It  shoulu  be  noted  that  c  is  not  exactly  the  speed  of 
sound.  Thus,  from  Eq.  (3) 


(dp/dp), j, 


=  nl 


(4) 


To  approximate  the  complete  UNESCO  function, 
we  have  determined  that 


p(S,  0,  p)  *  p(S,  0,  0)  +  I04  4  1  -  0.20  (5a) 

c-  cl 
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Modeling  Deep  Convection  in  the  Greenland  Sea 

S.  Hakkinen,'  G.  L.  Mellor.  and  L.  H.  Kantiia: 


Atmospheric  and  Oceanic  Sciences  Program.  Princeton  University ,  Princeton.  A V *  Jersey 

The  development  of  deep  convective  events  in  the  high-latnude  ocean  is  studied  using  a  three, 
dimensional,  coupled  ice-ocean  model.  Oceanic  mixing  ts  described  according  to  the  level  2.5 
turbulence  closure  scheme  in  which  convection  occurs  in  a  continuous  way,  i.e..  convecuve 
adjustment  is  not  invoked.  The  model  is  forced  by  strong  winds  and  surface  cooling.  Strong  upwciimg 
at  the  multiyear  ice  edge  and  consequent  entrainment  of  warm  Atlantic  waters  into  the  mixed  layer  is 
produced  by  winds  parallel  to  the  ice  edge.  Concomitant  cooling  drives  deep  convection  and  produces 
chimneylike  structures.  Inclusion  of  a  barolropic  mean  flow  over  topography  to  the  model  provides 
important  preconditioning  and  selects  the  location  of  deep  convection.  The  most  efficient  precondi¬ 
tioning  occurs  at  locations  where  the  flow  ascends  a  slope.  In  a  stratified  environment  similar  to  the 
Greenland  Sea  with  a  12  ms"1  wind  the  model  simulations  show  that  localized  deep  convection  takes 
place  after  about  10  days  to  depths  of  1000  m. 


I.  Introduction 

Deep  water  formation  at  high  latitudes  in  the  world  oceans 
provides  cold  water  masses  which  spread  equatorward  and 
balance  the  advection  of  the  warm  surface  waters  from  the 
low  latitudes.  Deep  water  formation  has  long  been  thought  to 
occur  on  the  continental  shelves  around  Antarctica  and  the 
Arctic  Ocean  as  a  result  of  the  brine-enriched  cold  shelf 
waters  sinking  to  the  deep  basins.  However,  it  has  not  been 
established  how  deep  and  on  what  time  scale  this  sinking 
occurs.  The  geostrophic  constraint  will  force  the  heavy 
waters  to  follow  the  coast,  and  only  the  frictional  forces  at 
the  bottom  can  cause  cross-isopycnal  (low.  Tracer  studies  in 
the  Arctic  show  that  the  deep  waters  in  the  Canadian  basin 
are  at  least  700  years  old  [Ost/tsnd  et  a!.,  1987],  which  means 
that  the  cold,  brine-enriched  shelf  waters  from  the  vast  East 
Siberian  shelves  have  not  been  sinking  to  the  deeper  basin. 
On  the  other  hand,  the  data  from  the  Eurasian  side  of  the 
Arctic  show  that  in  the  Nansen  basin,  the  deep  waters  are 
only  a  few  decades  old  [Wallace  and  Kryselt,  1987],  Osflttnd 
et  at.  [1987]  give  a  more  specific  age  of  30  years  for  the 
Eurasian  deep  water. 

Another  deep  water  source  is  open  ocean  production  via 
deep  convection,  most  notably  to  the  Greenland  and  Iceland 
seas  in  the  northern  hemisphere  and  to  the  Weddell  Sea  in 
the  southern  hemisphere.  In  the  southern  hemisphere  the 
manifestation  of  a  convection  event,  as  a  deep,  nearly 
uniformly  stratified  water  column  beneath  the  summer 
mixed  layer,  has  been  found  in  the  Weddell  Sea  by  Gordon 
(1978).  This  Weddell  Sea  "chimney’'  was  about  13  km  in 
diameter,  about  4-3  times  the  typical  baroclinic  Rossby 
radius  of  deformation  in  the  area.  There  is  only  one  known 
deep  convection  area  in  the  temperate  zone,  namely,  the 
Gulf  of  Lions  in  the  Mediterranean,  where  the  cold  northerly 
Mistral  and  Tramontaine  winds  cause  strong  mixing  and 
deep  convection  in  the  cyclonic  gyre  located  in  the  area 


'Now  at  NASA  Goddard  Space  Flight  Center.  Greenbelt.  Mary¬ 
land. 

'Now  at  Colorado  Center  for  Astrodynamics  Research.  Univer¬ 
sity  of  Colorado,  Boulder. 
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[ MEDOC  Group.  1970].  Deep  convection  extending  to  1500 
m  can  also  occur  in  the  Labrador  Sea  [Clarke  and  Gascard. 
1983]  where  this  process  is  responsible  for  the  formation  of 
the  Labrador  Sea  water,  an  important  component  of  the 
Atlantic  Jniermediate  water  [McCartney  and  Talley.  1984] 

Investigators  searching  for  chimneys  in  the  Greenland  Sea 
were  unable  to  find  these  mesoscale  features  unit!  very 
recent  years,  when  field  programs  such  as  the  Marginal  Ice 
Zone  Experiment  (MIZEX)  and  lhe  Greenland  Sea  projecl 
were  conducted.  During  lhe  winter  MIZEX,  some  convec¬ 
tive  areas,  about  10  km  across,  were  found  near  the  ice  edge 
region  [MIZEX  'S7  Group.  1989].  These  chimneys  extended 
at  least  to  250-300  m,  which  was  the  depth  limit  of  the 
(towed  Sea  Soar)  hydrographic  survey.  The  temperatures 
and  salinities  in  the  chimneys  ranged  from  0"  to  0.5*C  and 
from  34.84  to  34.94  ppt  respectively.  Rudels  et  at.  [1989] 
reported  finding  a  convective  area  wiih  a  depth  of  1230  m 
and  an  unknown  spatial  extent  from  a  1987  winter  cruise  to 
the  Greenland  Sea.  Bogorodsky  et  al.  [1987]  have  reported 
convective  areas,  about  100  km  across,  in  the  Greenland  Sea 
which  extended  to  the  bottom  depth  of  3500  m.  Johonrtessen 
et  al.  [1990]  have  found  further  evidence  of  wintertime 
convection  in  the  Boreas  basin,  the  northern  most  part  of  the 
Greenland  Sea.  In  all  of  the  aforementioned  studies  the 
information  of  the  horizontal  extent  of  the  chimneys  beyond 
the  observed  cross  section  is  uncertain.  Based  on  the  obser¬ 
vations  of  Johannessen  et  al.  [1990],  some  of  the  "chim¬ 
neys"  may  actually  be  elongated  features,  even  though  they 
suggest  that  chimneys  originate  from  cyclonic  eddies  which 
are  abundant  in  that  particular  area. 

On  the  basis  of  the  hydrographic  data  of  the  deep  convec¬ 
tive  area  by  Gordon  [1978],  Killworth  [1979]  used  a  nonpen- 
etrative  convection  model  to  study  possible  mechanisms  and 
necessary  surface  forcing  to  produce  these  features  which  he 
termed  "chimneys."  He  concluded  that  there  must  be  some 
preconditioning  mechanism  which  selected  a  particular  over¬ 
turning  area.  As  a  possible  mechanism  he  suggested  ba¬ 
roclinic  instability;  eddies  resulting  from  the  instability  occur 
at  the  same  spatial  scales,  and  the  pycnocline  in  the  cyclonic 
eddies  is  nearer  to  the  surface  than  the  surrounding  water, 
thus  giving  rise  to  enhanced  mixing. 

Martinson  et  al.  [1981  ]  later  revised  the  model  of  Kilhvorlh 
[1979]  to  study  deep  convection  in  (he  large  polynyas  occur- 
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Results  from  modeling  studies  of  the  ice-ocean  system  in  the  Arctic  Basin  and  in  the  Norsvegian- 
Greenland-Barents  seas  are  presented.  We  used  a  three-dimensional  coupled  ice-ocean  model 
developed  at  Princeton  University.  The  ocean  mode!  applies  the  primitive  equations  and  a  second 
moment  turbulence  closure  for  turbulent  mixing.  The  snow-ice  model  uses  a  three-level  thermody¬ 
namic  scheme  which  resembles  Scanner's  (1976a)  model.  Our  conclusions  based  on  the  seasonal 
simulations  are  as  follows.  (1)  Using  monthly  climatological  surface  heat  flux  and  wind  stress,  the 
seasonal  variability  of  the  ice  cover  is  quite  realistic  in  that  the  thickest  ice  is  located  north  of 
Greenland  and  the  average  ice  thickness  is  about  3  m.  The  largest  deviation  between  the  simulated  and 
observed  ice  cover  is  in  the  Greenland  Sea  where  oceanic  conditions  determine  the  ice  edge.  Basically, 
the  monthly  climatological  forcing  does  not  result  in  strong  enough  mixing  to  bring  sufficient  heat  from 
the  deep  ocean  to  keep  the  central  Greenland  Sea  gyre  ice  free.  The  results  improve  for  both  the  ice 
cover  and  ocean  by  invoking  daily  wind  forcing  for  which  we  first  chose  year  1987.  In  the  ocean  model, 
the  large  mixing  events  associated  with  storm  passages  are  resolved,  and  as  a  result,  the  overall 
oceanic  structure  in  the  Greenland  Sea  appears  to  be  more  realistic.  However,  no  deep  convection 
takes  place  in  the  model  during  1987  which  is  likely  the  result  of  diminished  storm  activity  in  the 
northern  part  of  the  Greenland  Sea.  The  ice  thickness  field  appears  to  be  very  anomalous  1987.  so  an 
experiment  with  1986  daily  wind  forcing  was  also  done,  which  resulted  in  an  ice  thickness  field  similar 
to  some  reported  from  other  ice  models.  (2)  Both  monthly  and  daily  surface  forcing  result  in  a  simdar 
behavior  of  the  Atlantic  waters  in  the  Arctic  Basin.  The  Atlantic  waters  circulate  at  about  the  observed 
level,  between  400  and  600  m.  The  survival  of  the  Atlantic  waters  in  the  basin  depends  strongly  on  the 
heat  loss  through  the  ice  cover,  and  it  appears  that  too  much  heat  is  lost  on  the  Eurasian  side  through 
the  ice  because  the  simulated  Atlantic  waters  are  too  cool  by  about  0.2-0.3’C.  (3)  For  the  monthly 
climatology  case,  a  large  amount  of  cold  and  salty  water  enters  the  Eurasia  Basin  from  the  Kara  and 
Laptev  seas  area  and  finds  its  way  toward  the  Canada  Basin.  This  water  mass  appears  to  result  from 
ice  formation  in  the  Kara  and  Laptev  seas.  When  applying  the  daily  forcing,  this  deep  salinity 
maximum  disappears  due  to  increased  mixing  on  the  shelves.  Nevertheless,  this  suggests  a  mechanism 
within  the  Arctic  Ocean  as  to  why  the  deep  Canada  Basin  is  much  saltier  than  the  Eurasia  Basin. 


1 .  Introduction 

The  objective  of  this  study  is  to  explore  the  circulation  of 
ice  and  water  masses  in  the  Arctic  Ocean  and  its  peripheral 
seas,  and  the  role  of  ice  thermodynamics-dynamics  in  these 
processes.  The  Arctic  Ocean,  together  with  Greenland, 
Iceland,  Norwegian,  and  Barents  seas,  comprise  a  unique 
system  where  a  large  portion  of  the  world  ocean  deep  waters 
are  formed.  Inflows  and  outflows  in  this  area  such  as  river 
runoff,  Atlantic  and  Pacific  inflow,  and  tee  export  determine 
the  overall  stratification  structure.  The  freshwater  compo¬ 
nent  is  especially  important  because  of  its  stabilizing  effect; 
in  the  Arctic  Ocean,  it  prevents  heat  exchange  between 
upper  and  deeper  parts  of  the  water  column;  in  the  Green¬ 
land  Sea,  an  excess  freshwater  cap  in  the  form  of  ice  can 
prohibit  the  renewal  of  deep  waters. 

The  Greenland,  Iceland,  Norwegian,  and  Barents  seas 
form  a  gateway  to  the  Arctic,  where  the  Atlantic  and  Arctic 
origin  waters  mix  and  form  new  water  masses  through 
convective  processes.  The  East  Greenland  Current  repre¬ 
sents  the  major  outflow  from  the  Arctic  Basin  and  it  is  a 
boundary  current  closely  following  the  continental  slope, 
The  Atlantic  waters  enter  the  area  through  the  Faeroe- 
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Shetland  Channel  and  flow  northward  along  the  Norwegian 
margin.  Part  of  this  current  branches  off  into  the  Barents 
Sea,  The  main  body  of  the  Norwegian  Current,  however, 
continues  northward  west  of  Svalbard,  where  it  mixes  with, 
and  plunges  beneath,  the  Arctic  surface  waters.  This  deep 
current  forms  a  counterclockwise  flow  in  the  Eurasia  Basin. 
The  same  water  mass  aiso  has  an  entrained  westward  branch 
which  forms  a  large  cyclonic  gyre  in  the  Greenland  Sea.  A 
schemetic  picture  of  the  different  surface  currents  is  shown 
in  Figure  1.  A  major  subsurface  outflow  from  the  GreenJand- 
Iceiand-Norwegian  seas  into  the  Atlantic  occurs  in  the 
Denmark  Strait  and  across  the  ridges  between  Iceland  and 
Scotland.  This  overflow  has  been  the  subject  of  intensive 
investigations  sponsored  by  the  International  Council  for  the 
Exploration  of  the  Sea  (ICES)  [Johannessen,  1986:  Swift, 
I986J. 

Compared  to  the  Greenland  Sea  our  knowledge  of  the 
Arctic  circulation  is  more  sketchy,  especially  for  the  subsur¬ 
face.  The  Arctic  surface  water  circulation  is  reasonably  well 
known  and  can  be  monitored  by  the  drift  patterns  of  sea  ice 
[Thorndike  and  Colony,  1982J.  The  surface  water  in  the 
Canada  Basin,  as  a  long-term  mean,  moves  in  an  anticy- 
clonic  gyre,  while  the  Eurasia  Basin  is  dominated  by  the 
Transpolar  Current  [Cordienko  and  Laktionov,  1969,  Figure 
1),  originating  in  the  Chukchi  Sea  and  flowing  out  in  Fram 
Strait.  The  circulation  of  the  mixed  layer,  however,  does  not 
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